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Nomenclature 

b = exposed electrodes width, 1mm 

c = clearance distance of consecutive electrodes, 0mm 

' c&f. = aerodynamic force coefficient in z-axis, F z /(qS) 

Ci p = roll moment coefficient caused by plasma actuators, 

d = chord of airfoil, 300mm 

e = insulated electrodes width, 4mm 

= clearance distance of consecutive plasma actuators, 6mm 

F, = aerodynamic force in z-axis, N 

, Qf r ,I y ,I z ) = the moment of inertia, kg ■ m 2 

(Ixy,Ixz, Iyz)= product of inertia, kg ■ m 2 

= total roll moment, N ■ m 

= aerodynamic roll moment, N ■ m 

= plasma-induced roll moment, N ■ m 

= coordinate origin 

= angular velocity in body axes, deg/s 

= dynamic pressure, ^pU^ 

= Reynolds number based on d and Uqc 

= area of the measuring section, mm 2 

= speed of freestream, m/s 

= wind span of airfoil, 3000mm 

= angle of attack, deg 

= air density, 1.225kg/m 3 

52</>, 0, ip) = Euler angles, deg 
Subscripts 

At = target value of A 
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1 Introduction 



PLASMA actuators, operating in atmospheric pressure air, have attracted increasing research interest in aerospace over 
the past two decades [U[2]. Most previous publications have focused on flow control applications that demonstrate the 
capability of plasma actuators and uncover, at least partially, the related fluid mechanics. In the present paper, we go one 
step further, studying the integration of flight control and flow control using plasma actuators, which constitutes the main 
contribution of this paper. As a demonstration, rolling maneuver of an airfoil was controlled only using plasma actuators, 
which could save mechanical moving parts of ailerons. In addition, an optimal flight controller was designed and demonstrated 
in simulations, taking account of flow control characteristics of plasma actuators. The proposed control method can also be 
considered for other flow control applications with plasma actuators. 

Various plasma actuators have been developed to address various flow control issues, which include high speed flow 
control using localized arc filament plasma actuators [3] and surface impulse discharges flat plate boundary layer flow 
control by non-thermal direct current (DC) corona [5] and aerodynamic-generated noise control with glow discharges [BJ. 
A comprehensive review of those plasma actuators can be found in the literature [2] and references therein. In this work, 
we adopted dielectric barrier discharge (DBD) plasma actuators [7] to control low speed aerodynamics and, going one step 
further, to control flight maneuver. A similar investigation has been conducted in a recent work [5], but using piezo-fluidic 
actuators. Satisfactory flight control performance has been shown in interesting flight tests jSJ. The present work differs 
from previous publications, focusing on the development of optimal flight control methodology for DBD plasma actuators. 



The principle of using DBD plasma actuators for flow control is not new. Figure 1(a) shows that a DBD plasma actuator 
normally consists of two electrodes isolated by a dielectric material. Potential candidates of the dielectric material include 
silicon rubber and flame retardant. An alternating current (AC) power supply (Fig. 1(b) I is applied to the two electrodes, 



between which glow discharges are generated, leading to the momentum transfer from charged nitrogen/oxygen particles to 
the local neutral gas through collisions. The collective effect is an induced fluid motion along the surface of the dielectric 
material, from the exposed electrode to the insulated electrode. The induced fluid motion develops vertical structures and 
manipulates local and global fluid mechanisms [5]. Detailed discussion of the related plasma and aerodynamic physics can 
be found in the literature [1]. 
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Figure 1: The plasma actuator, where (a) the schematic of a DBD plasma actuator; [TU] and (b) the AC voltage applied to 
the electrodes (top) and the plasma potential (bottom) measured in the plasma at 10 mm from the exposed electrode. 

In previous works, deployable flow effectors have been adopted on the upper surface of an NACA 0020 airfoil to delay flow 
separation [TT]. It has been shown that the stall angle can be increased from 18deg to 20deg. In addition to aerodynamic 
control, one natural extension is to control airfoil flight dynamics using plasma actuators. In particular, the lift modification 
due to plasma can generate desired pitching and rolling moments. As a result, it is possible to save the mechanical moving 
parts for aileron and flaps. However, feedback control system has to be considered along with the plasma active flow control 
to achieve optimal performance. A variable structure feedback control system has been proposed previously for cavity flow- 
induced noise control using glow discharges [12]. The fundamental idea behind that feedback control case was to adjust the 
authority of plasma actuation to a required strength, according to the feedback measurements and the variable structure 
model. The modulation of plasma actuation in real-time |13j . however, needs complex electrical circuits and generates serious 
electromagnetic pollutions. A simpler configuration of plasma actuators that only has two working states (on and off) was 
considered in this work. A bang-bang controller, also known as an on-off controller that abruptly changes actuation between 
upper and lower bounds, was used to feedback control the roll of NACA 0015 airfoil. The bang-bang control method is 
expected to generate less electromagnetic interference and is, inherently optimal in terms of its capability to achieve control 
objective in minimum time. 



2 



It can be seen from Fig. 1(b) that the AC input and the plasma potential have oscillations, which suggest the "on" 



state of a plasma actuator changes over time. It has been discovered in our previous work |10) that the maximum amplitude 
of the plasma-induced speed rapidly varies between 6m/s and 8m/s, leading to an oscillating lift force and force moment. 
The situation will be worsened for practical cases with time varied power supply at various atmospheric environments. This 
practical issue was considered and the disturbance rejection capability of the control method was studied in this work. 
Moreover, most DBD plasma actuators have a quite limited control performance for high speed flow cases. Different plasma 
actuators, such as impulse discharges [J] and corona [5] have been suggested for such cases. The attention of this paper 
focuses on low speed aerodynamic cases (around 15m/s) and can find applications in unmanned aerial vehicles. However, it 
is worthwhile to emphasize that the proposed control method is generic and could be applicable for applications with other 
plasma actuators. 

The following of this paper is organized as follows. Section II describes experimental apparatus that achieve aerodynamic 
data for the next flight control simulation. The bang-bang control method was briefly introduced in Sec. Ill to complete 
the paper. A bang-bang controller was thereafter designed particularly for the airfoil roll control case with DBD plasma 
actuators. The integration of flight control and active flow control was simulated in Sec. IV, where the oscillations of plasma 
actuations and the uncertainties in aerodynamic data were considered and addressed. A brief summary was provided in the 
end of this paper. 



2 Experimental Apparatus 



Roll control of an NACA 0015 airfoil (0.3m chord, 3m span) with plasma actuators was studied in this work. Figures [2] shows 
the body axes used throughout the rest of the paper. The origin O of the coordinates is located at the center of gravity of the 
airfoil. Plasma actuators are placed on the upper surface. The span of the plasma actuators is 0.75m, evenly covering the left 
and right wings (Fig. 0). Figure [5] shows an enlarged view of one plasma area, which is composed of nineteen consecutively 
spaced plasma actuators. The structure and configuration of each plasma actuator are similar to those in Fig. 1(a) All 
nineteen plasma actuators are parallel circuits and can be simultaneously switched on/off to a high voltage AC power supply. 
In roll control case, the plasma actuators on either left or right wing can be independently switched on/off. 



Plasma actuators 
3000 




Figure 2: The schematic of the NACA 0015 airfoil and the layout of plasma actuators for the roll control simulation case 
(dimensions are in millimeters). 

The aerodynamic data requested by flight control was achieved by conducting experiments in the FL-5 wind tunnel at 
Aerodynamics Research Institute. The FL-5 low-speed wind tunnel has an open-return circuit design. The open testing 
section is round with a diameter of 1.5m. The wind speed ranges from Om/s to 53m/s. Figure [3] shows the experimental 
setup. 

An NACA 0015 airfoil (0.3m chord, 0.9m span) was manufactured by Teflon and installed in the open test section. The 
rotating plates that can rotate within ±24deg are used to hold the model and test equipments as well as to maintain a good 
flow quality. Only the central section (0.2m span) of the model is covered with plasma actuators. As a result, the potential 
influence on the plasma flow control from the three-dimensional fluid and boundary flow local to the rotating plates can be 
omitted. The thickness of the dielectric material (epoxy polymer) is 1.5mm. The surface of the NACA 0015 airfoil can be 
etched to smoothly contain the plasma actuators. The electrodes are made by copper. The width of the exposed copper 
electrodes is b = 1mm. The width of the insulated electrodes is e = 4mm. The clearance distance c between two electrodes 
is zero. The clearance distance between two plasma actuators is / = 6mm. It should be pointed out that the above values 
are empirically chosen, reflecting a tradeoff between the experimental geometry limitation and plasma performance. It can 
be seen that the simulation setup shares the same geometrical setup for plasma actuators. Hence, the plasma induced roll 
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Figure 3: The setup in the wind tunnel with plasma glow discharges. 



momentum can be derived from the aerodynamic force measurements, although the experimental model is different from the 
simulation model for the convenience of the experiments. The force and moment coefficients of the simulation model, with 
and without plasma actuation, were respectively calculated based on experimental results, taking geometry difference into 
account. The AC voltage applied is 22kV, 180W at 4.7kHz. The following simulations adopt experimental results achieved 
at 15m/s. The corresponding Reynolds number with respect to the chord length is about 10 6 . 
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Figure 4: (a) Experimental results of C' z at Uoc = 15m/s, where the corresponding Reynolds number is 10 6 ; and (b) the 
related plasma-induced roll moment coefficient in the simulations. 



The experimental measurements of aerodynamic force are shown in Fig. 4(a) The force coefficient C' z = F z j (qS), where 
F z is the force along the z-axis, S is the area of the airfoil, q is dynamic pressure that is pU^/2, p is the air density, Uoo is 
the freestream speed. It can be seen that the amplitude of C z increases almost linearly along with the increase of a. The 
C z quickly drops beyond the so-called stall angle. Plasma actuation can slightly increase stall angle by 3deg. On the other 
hand, the increase of C z due to plasma actuation at low angles of attack is quite small. For example, we can only have 3.56% 
increase (from 0.705 to 0.73) in \C Z \ at a = 12deg. Figure [4(b) | shows the plasma-induced roll moment coefficient that is used 
in the simulations, where plasma actuators in the left airfoil are activated and the right ones are deactivated. The relatively 
large roll moment around 20deg is caused by the increased stall angle using plasma actuation. The roll moments at other 
angles are small. However, it can be seen in the following simulations that even such small roll moments can achieve good 
roll control performance. 
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3 Bang-Bang Controller Design 



In control theory, bang-bang control is well known as the minimum-time optimal feedback control method, whose control 
inputs are constrained to only two levels |14) . In particular, the roll dynamics of the airfoil can be described for the airfoil 
roll control case by 



L = I x p- (I y - I z )qr + I xy (pr - q) 
L = Lq + L p , 
p = <f> — (sin0)-0, 



I xz (pq + r) + I yz {r 2 - q 2 ) 



(1) 



where [I x , I y , I z ) is the moment of inertia in body axes, I xy is the product of inertia about ox and oy axes, I xz is the product 
of inertia about ox and oz axes, I yz is the product of inertia about oy and oz axes. The angular velocities are represented 
by (p,q,r). The Euler angles of the airfoil about the flat earth are (<j>,6,i/j). L is the total roll moment, which consists of 
aerodynamic roll moment Lq and plasma- induced roll moment L p . 

The produces of inertia I xy = I yz = 0, as the airfoil is symmetric about the oxz plane and the airfoil has a uniform mass 
distribution. In addition, the pitch and yaw motions are excluded for simplicity, i.e. q = 0, r = 0, ip = 0, 9 = 0. As a result, 
4> = p, where p is the roll angular velocity in body axes. 

It is implicitly assumed that the airfoil is already trimmed without plasma actuation. The plasma-induced moment is 
L p = M in case the plasma actuators on the left wing are activated. Similarly, L p — —M if the plasma actuators on the right 
wing are activated. The scenario that all plasma actuators on both wings are activated was not considered in this work. In 
summary, the roll dynamics can be simplified to 

$ = ±M/I X . (2) 
A state space model can be accordingly formulated as 
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(3) 



where the control input u is the sign of plasma- induced roll moment. 

With no loss of generality, the initial state of the airfoil is assumed at [</>(0), 0(0)] = [0, 0]. The objective of the roll control 
is [4>(t), 4>(t)] = [<pT, 0], where 4>t is the target roll angle. It is straightforward to achieve 
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for u = 1: I fir for u = — 1, 



(4) 



where C\ is a constant value. Figure 5(a) shows the phase line (f> = ± 2 a//j (0) 2 + = which is achieved by 

eliminating the variable t in Eqs. (H]). 




Figure 5: Schematics of bang-bang control, where (a) phase line with u = ±1; (b) the on-off line for roll control. 



As shown in Fig. 5(b)| the initial state [0(0), 0(0)] = [0,0] is at point O. The control objective [0(f), 4>(t)} = [0t,O] at B. 
As a result, the initial control input should be set to u = +1 to turn the initial state O along the phase line <j> = 2 m/i 2 ^° 

state A, which is the intersection on the other phase line <f> = — 2 m/i T (0 ■ > 0). The control input switches to u = — 1 
once A is reached. The state thereafter advances to the target state B[0t,O]. 
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Figure 6: The simulation case in MATLAB. 
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Figure 7: Dynamics during the airfoil rolls from Odeg to lOdeg at a = 12deg and Uoo — 15m/s, where (a) roll angle dynamics 
and (b) phase line. 



A so-called on-off line (the dashed line in Fig. 5(b) I only depends on <f>T and can be formulated as 

1 



In summary, for any initial solution, the desired control input it is: 

1 . 1 



(5) 
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(6) 



4 Simulation and Discussion 

Numerical simulations were conducted to demonstrate airfoil roll control using plasma actuators. Figure [6] shows the imple- 
mentation of the simulation case in MATLAB. The block of bang-bang controller implements Eqs. (H2). The wing dynamic 
block implements the airfoil roll dynamics. The moment of inertia I x is 4.05kgm 2 , given the span length (3m) and the 
uniformly distributed mass (5.4kg) with the density of 200kg/m 3 . The dynamic pressure is 137. 8Pa, given the air density 
1.225kg/m 3 and mean flow velocity 15m/s. 

Figure [7] shows the simulation results for <f>T = lOdeg, where the angle of attack a = 12deg, at which it is implicitly 
assumed that the airfoil is already trimmed. It can be seen that plasma-induced moments roll the airfoil to the target angle 
in about 1.5s. Almost no overshoot can be found. The satisfactory control results with the bang-bang control method 
may be arguable as the plasma- induced aerodynamic inputs are varying over time (as suggested by Fig. 1(b)). Figure [5] 
is the plasma-induced velocity measured at 20mm downstream from the exposed electrodes [5]. It can be seen that the 
plasma-induced velocities oscillate at high voltage input (Vpp=22kV and 12.5kV). The variance is approximately 20%. 
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Figure 8: Measurements of plasma-induced velocities in the stationary atmospheric pressure air. 
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Figure 9: Monte Carlo simulation results of roll control from Odeg to lOdeg, where the plasma-induced roll momentum is 
randomly chosen, where (a) roll angle dynamics and (b) phase line. 



A Monte Carlo simulation was conducted to examine the effect of the varying plasma actuations on the roll control. The 
plasma-induced momentum is approximated with ±20% uncertainty about the nominal values Ci p in Fig. |4(b)| Figures IMTUl 
show the Monte Carlo simulation results. The roll moment coefficient is randomly chosen between 0.8Cz p and 1.2C; p in 100 
repeated simulations. It can be seen in Fig. [5] that overshoots could appear for varied plasma-induced momentum. The 
phase line slightly changes as well. In addition, the rising time of the control has been affected. 

Figure [TU1 quantitatively shows the changes. The roll-mode time constants (the time needed to reach 63.2% of the target 
roll angle) for different roll targets (from lOdeg to 90deg) at a 



12deg are shown in Fig. 10(a) 



According to MIL-F-8785C 

specification of flying quality, the maximum roll-mode time constant is 10s for level 3 flights of light airplanes. Figure [10] 
shows that this flying quality section can be satisfied with the proposed bang-bang controller and plasma actuators. Figure 
|10(b)| shows the roll-mode time constants for 10 roll degree at various angles of attack between 4deg and 24deg. It can be 
seen that the roll control is most effective at a high angle of attack around 20deg, where the flow separation has been delayed 
with plasma actuation. The flying qualities (< 10s) are still satisfied. 

It is more practical to regard the almost 20% aerodynamic variance as an instantaneous disturbance to the controlled 
system. The Fourier spectrum of the results in Fig. [8] largely lies between 1Hz and 30Hz. As a result, the plasma-induced 
roll moments are varied at these two frequencies (Fig. ll(a)|ll(d) I to verify the robustness of the proposed control method. 

show that slight change can be found in roll dynamics and phase line (< 5%) for the 1Hz perturbation 



11(c) 



11(e) 



11(f) suggests the control method is insensitive to the relatively high frequency disturbance at 30Hz. 



Figures 11(b) 
case. Figure 

Almost no change can be found in the roll dynamics and phase line. As a result, the proposed bang-bang control is robust 
and can address practical issues of plasma actuators. 

At last, Fig. Q2] shows the roll dynamics of the airfoil commanded by a series of roll commands. The nominal roll angles 
with constant plasma-induced roll moments satisfactorily track the roll commands with less than 0.5s time delay. On the 
other hand, the perturbed roll angles with 20% variances of moments at 1Hz also follow the roll commands well, suggesting 
the good performance of the proposed control method. 
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Figure 10: Roll-mode time constant error bar where (a) a — 12deg; (b) (f> = lOdeg. 




Figure 11: Simulation results of roll control from Odeg to lOdeg at a = 12deg, where (a)(d) 20% variance of plasma- induced 
roll moments at 1Hz and 30Hz, (b)(e) the corresponding roll angle dynamics, and (c)(f) the related phase lines. 
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Figure 12: Roll angle dynamics for a series of roll commands. 



5 Summary 

The main contribution of this Note is to integrate flight control with active flow control using plasma actuators. The bang- 
bang control method has been proposed for plasma actuators, taking account of practical issues, such as limited actuation 
states with instantaneously varied aerodynamic control performance. Flow control effects have been examined in wind tunnel 
experiments, which show that the plasma authority for flow control is limited. Flow control effects are only obvious at pitch 
angles near stall. However, flight control simulations suggest that, using the proposed optimal control method, even those 
small plasma-induced roll moments can satisfactorily fulfill the maneuver tasks and meet flight quality specifications. In 
addition, the disturbance from volatile plasma-induced roll moments can be rejected well. Hence, the proposed bang-bang 
control method is a promising candidate of control design methodology for plasma actuators. The ongoing and future works 
include airfoil flight control in wind tunnel and final flight tests. 
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